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Minimalizacja efektu katalizy w metodach doboru zmiennych 


Minimalization of the Catalysis Effect in the Methods of the Choice of Variables 


Zagadnienie doboru zmiennych w jednoröwnaniowych liniowych mo- 
delach ekonometrycznych bylo przedmiotem uwagi licznych autorów. 
Stad w pismiennictwie naukowym istnieje wiele artykulów, jak tez opra- 
cowan ksigzkowych, z ktörych najwazniejsze to [4], [5], [11]. Do metod 
doboru zmiennych najczęściej stosowanych należą między innymi: meto- 
da selekcji a posteriori, metoda selekcji a priori, metoda regresji kroko- 
wej, metoda Allena, metoda Hellwiga (por. [4], metoda suboptymalnej 
wartości współczynnika integralnej pojemności informacyjnej (por. [9]), 
metoda Rogowskiego (por. [10]), metoda Guzika ([6]). Dobierając zmienne 
do liniowego jednorównaniowego modelu ekonometrycznego powinno mieć 
Się na uwadze takie postulaty jak: 

1) maksymalizację stopnia dokładności, z jaką model opisuje dane 
zjawisko, 

2) istotność parametrów strukturalnych modelu, 

3) dokładność ocen parametrów strukturalnych modelu, 

4) koincydencję modelu, 

5) brak zmiennych współliniowych, 

6) brak efektu katalizy lub minimalne jego natężenie. 

Metoda selekcji a posteriori, metoda selekcji a priori, metoda regresji 
krokowej, metoda Allena, metoda Hellwiga oraz metoda Guzika nie za- 
pewniają koincydencji modelu. W pracy [2] przeprowadzono adaptację 
czterech pierwszych metod w celu uzyskania przy ich pomocy modeli 
koincydentnych. Adaptacja tych metod polegała na wkomponowaniu do 
alogrytmów właściwych tym metodom warunku zapewniającym koincy- 
dentność modelu. Uzyskano to poprzez zastosowanie współczynników ko- 
relacji częstkowej do badania koincydentności, gdyż znak oszacowania pa- 
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rametru strukturalnego modelu jest identyczny ze znakiem współczynnika 
korelacji cząstkowej. 

Model otrzymany metodą Rogowskiego oraz metodą suboptymalnej 
wartości współczynnika integralnej pojemności informacyjnej jest koin- 
cydentalny. Żadna z wyżej wymienionych metod nie uwzględnia proble- 
mu efektu katalizy. 

W niniejszej pracy przedstawimy dwie metody doboru zmiennych, w 
wyniku których otrzymany model jest koincydentny o minimalnym na- 
tężeniu efektu katalizy. Efekt katalizy został wykryty i opisany przez 
Z. Hellwiga w pracy [7]. Poniżej przypomnimy definicję efektu katalizy. 


DEFINICJA 1 (por. [7]) 


Mówimy, że w modelu ekonometrycznym określonym przez regularną 
parę korelacyjną (R, Ro) występuje efekt katalizy, jeżeli istnieje taka 
para wskaźników (i, j), dla której 


r 
Ty <0 lub Ty > a 


gdzie R=[ry]k«k, Ro= [ri]kxı są odpowiednio macierzą korelacji dla zmien- 
nych objaśniających i wektorem, którego i-tą współrzędną jest współczyn- 
nik korelacji między i-tą zmienną objaśniającą a zmienną objaśnianą. Po- 
nieważ efekt katalizy może występować w modelu z różnym natężeniem, 
dlatego został określony miernik (por. [7]). 


n=r?—H, 
gdzie r? jest kwadratem wspólczynnika korelacji wielorakiej, zas H jest 
pojemnością integralną informacji służący do mierzenia natężenia efektu 
katalizy. Na podstawie nierówności 

0<n<1 
wykazanej w [8] widzimy, ze w modelu nie wystepuje efekt katalizy gdy 
n=0, natomiast natężenie efektu katalizy osiąga wartość największą, gdy 
n=1. W metodach, które zaproponujemy w niniejszej pracy wykorzysta- 
my następującą definicję 2 oraz twierdzenie 1. 

DEFINICJA 2 (por. [3]). 

Funkcję f określoną na parze korelacyjnej (R, Ro) i przyjmującą war- 
tość w przedziale <0,1>, która jest ciągła względem wszystkich elemen- 
tów macierzy R oraz wszystkich współrzędnych wektora R. oraz taką, że 

A k 2 = 
gdy R jest macierzą jednostkową stopnia k oraz Dizi =1, to f(R, Ro)=1 
nazywamy miernikiem liniowego jednorównaniowego modelu ekonome- 
trycznego. 
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TWIERDZENIE 1 (por. [3]) 


Jeżeli f, oraz f} sa miernikami modelu w snsie definicji 1, to dla do- 
wolnych liczb rzeczywistych dodatnich c, oraz cz funkeja f = fi fs? jest 
miernikiem modelu. 

Zauwazmy, ze w sensie definicji 1 miernikami modelu sa: 

— integralny współczynnik pojemności informacji H, 

— integralny współczynnik koincydencji M?=r? det R zdefiniowany 
po raz peirwszy w [1], 

— współczynnik K=1 —n=q2+H (9*=1 — r? jest współczynniem in- 
determinacji) mierzący natężenie efektu katalizy w ten sposób, że gdy 
w modelu brak jest efektu katalizy to K=1, natomiast wraz ze wzrostem 
natężenia katalizy współczynnik K maleje od 1 do 0; wartość 0 przyjmuje 
dla największego natężenia efektu katalizy. 

W myśl twierdzenia 1 z miernikow H, M? i K tworzymy następujące 
dwa mierniki: 

Hk=H'K, 
N=M? K 
które są miernikami modelu. Wiadomo, że miernik H przy doborze zmien- 
nych do modelu uwzględnia zarówno korelację między zmienną objaśnia- 
ną a zmiennymi objaśniającymi, jak też wzajemną korelację między 
zmiennymi objaśniającymi. Współczynnik H=1 wtedy i tylko wtedy, gdy 
r?=1 oraz R=I, gdzie I jest macierzą jednostkową. Wartość H=0 ozna- 
cza brak jakiejkolwiek korelacji między zmienną objaśnianą a zmiennymi 
objaśniającymi. Ponieważ w modelu powinny występować takie zmienne 
objaśniające, dla których wartości H, K i M? są możliwie duże dlatego 
pożądana jest również duża wartość miernika Hy, oraz N. 
Poniżej przedstawimy pierwszą z zapowiedzianych metod. Niech 


y= {X;, X2,..., Xx) będzie zbiorem potencjalnych zmiennych objaśnia- 
jących. Przez y; (i=1,2,..., 2*—1) oznaczmy i-ty podzbiór zbioru x. Dla 
każdego podzbioru y (i=1,2,..., 2*—1) obliczamy współczynnik Hx 


oznaczając go symbolem Hy (xi). Następnie postępujemy według algo- 
rytmu: 


1) Znajdujemy zbiór xi, spełniający warunek 
Hi(xu) = ax (Hr ba) eC KX) 
2) Sprawazamy, czy model o zmiennych objasniajacych nalezacych 


do zbioru Xi: jest koincydentny. Jesli tak to procedura jest zakonczona i do 
modelu zostają wybrane zmienne objaśniające ze zbioruxi. Jeśli nie to 
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3) Wyznaczamy zbiór xi,taki, że 


Hy (Xi2) = ‚max He) xi AK In COX) 


4) Sprawdzamy, czy model, w ktörym zmienne objasniajace naleza do 
zbioru xi, jest koincydentny. Jesli tak, to konczymy algorytm i za zmien- 
ne objasniajace do modelu przyjmujemy zmienne ze zbioruxi. Jeśli nie, 
to 


5) Wyznaczamy zbiór xisspelniajacy warunek 
Ax(xis) = max (Hi): xi A Xn IM A Xi i xi Cx] 
1<i52k--3 


6) Postepujemy analogicznie jak w 4) itd. 
Koincydentnost modelu sprawdzamy wykorzystujac nastepujacy wa- 


runek konieczny i wystarczajacy na to, by dowolna zmienna objasniajaca 
miała własność koincydencji. 


TWIERDZENIE 2 (POR. [9]) 


Zmienna objaśniająca X, jest koincydentna w modelu określonym 

przez regularną parę korelacyjną (R, Ro) wtedy i tylko wtedy, gdy 
r > [Ri] [Ru] "R 

gdzie R; oznacza i-tą kolumnę macierzy R z pominięciem i-tej wspölrzed- 
nej, R, powstaje z wektora R» przez odrzucenie i-tej współrzędnej, na- 
tomiast Ri oznacza podmacierz otrzymaną z macierzy R przez skreślenie 
wiersza oraz kolumny o numerze i-tym. 

Model jest koincydentny, jeżeli każda zmienna objaśniająca tego mo- 


delu ma własność koincydencji. W celu łatwiejszego korzystania z twier- 
dzenia 2 proponujemy wykorzystać następujące. 


TWIERDZENIE 3 


Jeżeli na elementach macierzy A; postaci 


SCE 


gdzie A;=[ay], t= [f], g=[gi] (,j=1,2,....,k), ZER (R jest zbiorem 
liczb rzeczywistvch) wykonamy takie przekształcenia elementarne (roz- 
poczynając od elementu a;,), które sprowadzają ją do postaci górnej ma- 
cierzy trójkątnej D= [d;;], to zachodzi równość 


Z—gTrAf=dke, k+1. 
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dla macierzy A; postaci 


Ry Ri 
A; = | [Ri]? ri | 


Na elementach macierzy A; zgodnie z twierdzeniem 3 wykonujemy od- 
powiednie przekształcenia elementarne spowadzające ją do postaci ma- 


cierzy trójkątnej D;= lat] (P, g=1,2,...,k). Jeżeli dl) >0 to zmien- 
na X; jest koincydentna w modelu określonym przez parę korelacyjną 
(R, Ro). Jeżeli natomiast di) <0to zmienna X; nie jest koincydentna. Mo- 
del, ktöry otrzymamy w wyniku opisanej wyZej procedury ma najwiek- 
szą wartość miernika H, wśród modeli koincydentnych o zmiennych objas- 
niających należących do podzbiorów  (i=1,2,...,24—1). 

Druga proponowana metoda wykorzystuje integralny miernik koin- 
cydencji M? oraz następujące 


TWIERDZENIE 4 (POR. [1]) 


Zmienna objaśniająca X;(i=1,2,...,,k) w modelu określonym przez 
regularną pare korelacyjną (R, Ro) jest koincydentna, jeśli spełniona jest 
nierówność y 2 

M? > M; 


gdzie M; : jest integralnym miernikiem koincydencji modelu, który pow- 
staje z wyjściowego modelu przez odrzucenie j-tej zmiennej objaśniają- 


cej. W metodzie tej najpierw wyznaczamy wszystkie podzbiory zbioru 
X= {X X»,..., Ex} (jest ich 2*—1) i dla każdego z nich obliczamy war- 
tość miernika N. Wartość miernika N dla i-tego podzbioru y; oznaczamy 
symbolem N (y;). 

Procedura postepowania jest nastepujaca: 

1) spośród wartości N (xi) (i=1,2..., 28-1) wybieramy wartość naj- 
większą. Podzbiör zbioru x, któremu odpowiada ta największa wartość 
oznaczamy symbolem x; czyli mamy nierówność 


N(Xi1) 2 Nix), i= 12 242° -1 
2) Sprawdzamy, czy zachodzi nierówność 
A M? (xi) > M} (Xi), gdzie M? (xi) 
jest integralnym miernikiem koincydencji 
dla modelu o zmiennych objasniejacych nalezacych do zbioru xi, nato- 


miast M?(xi,) jest integralnym miernikiem koincydencji obliczonym dla 
modelu o zmiennych objaśniających ze zbioru xi, bez j-tej zmiennej ob- 
jaśniającej. Jeżeli tak, to postępowanie kończymy i za optymalny zbiór 
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zmiennych objasniajacych uznajemy zbiör xi, Otrzymany model ma wlas- 
nos¢ koincydencji w mysl twierdzenia 4. Jezeli nie, to 

3) Spośród wartości N (xi) dla i=1,2,...,2*—2 oraz i Æi; wybieramy 
oznaczamy symbolem x;, 
wartość największą. Podzbiór, któremu odpowiada ta en wartość 

4) Sprawdzamy, czy zachodzi nierówność 

AM?) > M} (Xi) 

(xi). i 
Jeżeli tak, to za optymalny zbiór zmıennych objasniających przy jmu- 
jemy zbiór xi, Jeśli nie, to znów wybieramy podzbiór xi, itd. 

Model, który otrzymamy w wyniku opisanej metody jest modelem 
o największej wartości miernika N spośród modeli koincydentnych. 

Poniżej przedstawimy przykład ilustrujący drugą z zaproponowanych 
metod. 

Przykład 
Niech y= {X,, X, X3). Zalözmy, ze dla zmiennych objasniajacych i zmien- 
nej objaśnianej macierze R i R, mają postać: 


1 03 05 0,6 
R=| 031 0,8 R=| 02 
0,5 0,8 1 0,3 


Zbiör x ma 7 nastepujacych podzbiorów: 


=(K, Xz, Xo}, w= (Ki, Xe}, w= (XX), = (Kaa KH}, GEO >. 

= (X), X {X3}. 
Bi każdego podzbioru zbioru x obliczamy wartość miernika N i otrzy- 
mu jemy: 
N (x1)= 0,08, N (x2)=0,30, N (x3)=0,25, N (x)= 0,04, 
N (x5) 0,36, N (4s)= 0,04, N (x)= 0,09. 

Spośród wartości N (xi) (i=1,2,...,7) wybieramy wartość najwięk- 
szą czyli 


N DSN (ys) 5 N(u)  G=1,2,...,7). 
Ponieważ podzbiór ys ma tylko jedną zmienną objaśniającą X,, więc mo- 
del opisywany przez tą zmienną jest koincydentny (model o jednej zmien- 
nej objaśniającej zawsze jest koincydentny por. ([2]). Stąd za otpymalny 
zbiór zmiennych objaśniających przyjmujemy zbiór ys. 
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SUMMARY 


The paper presents a proposition of such two methods of the choice of varia- 
bles to the linear one-equation econometric model which ensure the model coin- 
cidence, a minimum intensification of the catalysis effect, a strong correlation 
between the explanatory variables and the explained variable as well as a poor 
correlation between the explanatory variables. 


